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We propose a theory to explain the formation of spiral arms and of all types 
of outer rings in barred galaxies, extending and applying the technique used 
in celestial mechanics to compute transfer orbits. Thus, our theory is based on 
the chaotic orbital motion driven by the invariant manifolds associated to the 
periodic orbits around the hyperbolic equilibrium points. In particular, spiral 
arms and outer rings are related to the presence of heteroclinic or homoclinic 
orbits. Thus, i?i rings are associated to the presence of heteroclinic orbits, 
while R1R2 rings are associated to the presence of homoclinic orbits. Spiral 
arms and R2 rings, however, appear when there exist neither heteroclinic 
nor homoclinic orbits. We examine the parameter space of three realistic, 
yet simple, barred galaxy models and discuss the formation of the different 
morphologies according to the properties of the galaxy model. The different 
morphologies arise from differences in the dynamical parameters of the galaxy. 

1 Introduction 

Bars are a very common feature of disc galaxies. In a sample of 186 spirals, 56% 
of the galaxies in the near infrared are strongly barred, while an additional 
16% are weakly barred [9]. A large fraction of barred galaxies show two clearly 
defined spiral arms [8], departing from the end of the bar. This is the case for 
instance in NGC 1300, NGC 1365 and NGC 7552. Spiral arms are believed to 
be density waves in a disc galaxy [13]. In [22], Toomre found that the spiral 
waves propagate towards the principal Lindblad resonances of the galaxy, 
where they damp down, and thus concludes that long-lived spirals need some 
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replenishment. Danby argued that orbits in the gravitational potential of a 
bar play an important role in the formation of arms [6] and Kaufmann & 
Contopoulos argue that, in self-consistent models for three real barred spiral 
galaxies, spiral arms are supported also by chaotic orbits [11]. 

The origin of rings has been studied by Schwarz who calculated the re- 
sponse of a gaseous disc galaxy to a bar perturbation [18, 19, 20]. He pro- 
posed that ring-like patterns are associated to the principal orbital resonances, 
namely ILR (Inner Lindblad Resonance), CR (Corotation Resonance), and 
OLR (Outer Lindblad Resonance). There are different types of outer rings. 
Buta classified them according to the relative orientation of the ring and bar 
major axes [4]. If these two are perpendicular, the outer ring is classified as 
Ri . If they are parallel, the outer ring is classified as R2 ■ Finally, if both types 
of rings are present in the galaxy, the outer ring is classified as R\R2- 

In Romero-Gomez et al. [16, 17], we note that spiral arms and rings em- 
anate from the ends of the bar and we propose that rings and spiral arms are 
the result of the orbital motion driven by the invariant manifolds associated 
to the Lyapunov periodic orbits around the unstable equilibrium points. In 
Romero-Gomez et al. [16], we fix a barred galaxy potential and we study the 
dynamics around the unstable equilibrium points. We give a detailed defini- 
tion of the invariant manifolds associated to a Lyapunov periodic orbit. For 
the model considered, the invariant manifolds delineate well the loci of an rRi 
ring structure, i.e. a structure with an inner ring (r) and an outer ring of the 
type Ri. In Romero-Gomez et al. [17], we construct families of models based 
on simple, yet realistic, barred galaxy potentials. In each family, we vary one 
of the free parameters of the potential and keep the remaining fixed. For each 
model, we numerically compute the orbital structure associated to the invari- 
ant manifolds. In this way, we are able to study the influence of each model 
parameter on the global morphologies delineated by the invariant manifolds. 

Voglis, Stavropoulos & Kalapotharakos study the chaotic motion present 
in self-consistent models of both rotating and non-rotating galaxies, conclud- 
ing that rotating models are characterised by larger fractions of mass in chaotic 
motion [23]. Patsis argues that the spiral arms of NGC 4314 are due to chaotic 
orbits and, to show it, he computes families of orbits with initial conditions 
near the unstable equilibrium points [15]. Voglis, Tsoutsis & Efthymiopoulos 
[24] reproduce a spiral pattern found in a self-consistent simulation using the 
apocentric invariant manifolds of the short-period family of unstable periodic 
orbits. They give the angular position of the apocentres, which is where they 
state the stars spend a large part of their radial period, as a soliton-type 
solution of the Sine-Gordon equation. 

In Sect. 2, we first present the galactic models used in the computations 
and the equations of motion. In Sect. 3, we give a brief description of the 
dynamics around the unstable equilibrium points and the role the invariant 
manifolds play in the transfer of matter around the galaxy. In Sect. 4, we show 
the different morphologies that result from the computations. 



Formation of spiral arms and rings in barred galaxies 



3 



2 Description of the model and equations of motion 

In this section, we first describe the bar models used in the computations by 
giving the density distributions, or the potentials used. We then write the 
equation of motion and we define the efi^ective potential and Jacobi constant. 



2.1 Description of the model 

We use three difl[erent models, all three consisting of the superposition of 

an axisymmctric component and another bar-like. Our first model is that of 
Athanassoula [2]. The axisymmetric component is composed of a disc, mod- 
elled as a Kuzmin-Toomre disc [12, 21] of surface density a{r): 

yi / 2 \ -3/2 
and a spheroid modelled by a density distribution of the form p{j): 

/ 2 \ -3/2 

p{r) = p, (^1 + y • (2) 

The parameters Vd and set the scales of the disc velocities and radii, re- 
spectively, and ph and determine the concentration and scale-length of the 
spheroid. 

Our bar potential is described by a Ferrers ellipsoid [10] whose density 
distribution is: 

/ po(l - m2)" m < 1 

\0 m>l, ^"^^ 

where w? = jc? + y'^ /b"^. The values of a and h determine the shape of the 
bar, a being the length of the semi-major axis, which is placed along the x 
coordinate axis, and h being the length of the semi-minor axis. The parameter 
n measures the degree of concentration of the bar and po represents the bar 
central density. 

Wc also use two further ad-hoc bar potentials, namely a Dehnen's bar 
type, ^1, (Dehnen [7]): 



r < a 



<l^i{r,e) = -\evlcos{2e)l , (4) 



r > a. 



and a Barbanis-Woltjer (BW) bar type, ^2, (Barbanis & Woltjer [3]): 

<Z>2(r, 9) = e^f?{ri - r) cos(26') (5) 

The parameter a is a characteristic length scale of the Dehnen's type bar 
potential, and vq is a characteristic circular velocity. The parameter e is related 

to the bar strength. The parameter ri is a characteristic scale length of the 
BW bar potential and the parameter e is related to the bar strength. 
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2.2 Equations of motion 

In order to compute the equations of motion, we take into account that the 
bar component rotates anti-clockwise with angular velocity Jlp = X2pZ, where 
Op is a constant pattern speed The equations of motion in this potential 
in a frame rotating with angular speed ftp in vector form are 

j: = -V* - 2(f2p X f) - rip X (rip X r), (6) 

where the terms — 2np x r and —i^p x (ftp x r) represent the Coriolis and 
the centrifugal forces, respectively, and r is the position vector. Defining an 
effective potential: 

^,s = ^-\n'p{x' + y% (7) 

Eq. (6) becomes = — V$cfF — 2(Jlp x f), and the Jacobi constant is Ej = 
1 1 f 1^ + <?eff) which, being constant in time, can be considered as the energy 
in the rotating frame. 

3 Dynamics eiround Li and L2 

For our calculations we place ourselves in a frame of reference corotating with 
the bar and place the bar major axis along the x axis. In this rotating frame we 
have five equilibrium points, which, due to the similarity with the Restricted 
Three Body Problem, are called Lagrangian points. The points located on the 
origin of coordinates, namely L3, and along the y axis, namely L4 and L5, 
are linearly stable. The ones located symmetrically along the x axis, namely 
Li and L2, are linearly unstable. Around the equilibrium points there exist 
families of periodic orbits, e.g. around the central equilibrium point the well- 
known xi family of periodic orbits [5] that is responsible for the bar structure. 

The dynamics around the unstable equilibrium points is described in detail 
in [16]; here wc give only a brief summary. Around each unstable equilibrium 
point there also exists a family of periodic orbits, known as the family of 
Lyapunov orbits [14]. For a given energy level, two stable and two unstable sets 
of asymptotic orbits emanate from the periodic orbit, known as the stable and 
the unstable invariant manifolds, respectively. We denote by W*. the stable 
invariant manifold associated to the periodic orbit 7 around the Lagrangian 
point Li, i = 1,2. The stable invariant manifold is the set of orbits that tends 
to the periodic orbit asymptotically. Similarly, we denote by W^. the unstable 
invariant manifold associated to the periodic orbit 7 around the Lagrangian 
point Li, i = 1,2. The unstable invariant manifold is the set of orbits that 
departs asymptotically from the periodic orbit (i.e. orbits that tend to the 
Lyapunov orbits when the time tends to minus infinity) (Fig. 1). Since the 
invariant manifolds extend well beyond the neighbourhood of the equilibrium 
points, they can be responsible for global structures. 

Bold letters denote vector notation. The vector z is a unit vector. 
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Fig. 1. Unstable (in dark grey), W^^, and stable (in light grey), W^^ , invariant 
manifolds associated to the periodic orbit around Li (in white solid line). In black 
solid lines, we plot the zero velocity curves for this energy level and the dashed curve 
shows the outline of the bar. 

In [17], we give a detailed description of the role invariant manifolds play in 
global structures and, in particular, in the transfer of matter. Simply speak- 
ing, the transfer of matter is characterised by the presence of homoclinic, 
heteroclinic, and transit orbits. Homoclinic orbits correspond to asymptotic 
trajectories tp such that S PF". n W^. , i = 1,2. That is, they arc asymptotic 
orbits that depart from the unstable Lyapunov periodic orbit 7 around Li 
and return asymptotically to it (Fig. 2a). Heteroclinic orbits are asymptotic 
trajectories tp' such that tp' e W". fl W^., i,j = 1,2, i ^ j. That is, they 
are asymptotic orbits that depart from the periodic orbit 7 around Li and 
asymptotically approach the corresponding Lyapunov periodic orbit with the 
same energy around the Lagrangian point at the opposite end of the bar Lj, 
i ^ j (Fig. 2b). We are interested in the homoclinic and heteroclinic orbits 
corresponding to the first intersection of the invariant manifolds with an ap- 
propriate surface of section. There exist also trajectories that spiral out from 
the region of the unstable periodic orbit and we refer to them as transit orbits 
(Fig. 2c). These three types of orbits are chaotic orbits since they fill part of 
the chaotic sea when we plot the Poincare surface of section (e.g. the section 
{x, x) near Li). 

4 Results 

One of our goals is to check the influence of each main free parameter of the 
models introduced in Sect. 2. In order to do so, we make families of models in 
which only one of the free parameters is varied, while the others are kept fixed. 
Our results in [17] show that only the bar pattern speed and the bar strength 
have a considerable influence on the shape of the invariant manifolds and, 
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Fig. 2. Homoclinic (a), hotcroclinic (b) and transit (c) orbits (black thick linos) in 
the configuration space. In dark grey lines, we plot the unstable invariant manifolds 
associated to the periodic orbits, while in light grey we plot the corresponding stable 
invariant manifolds. In dashed lines, we give the outline of the bar and, in (b) and 
(c), we plot the zero velocity curves in dot-dashed lines. 



thus, on the morphology of the galaxy. Having established this, we perform 
a two-dimensional parameter study for each bar potential and we obtain all 
types of rings and spiral arms. 

In Fig. 3 we show the model rings and the spiral structure we obtain with 
our models. We plot the unstable (Fig. 3a, b and d) and the unstable and stable 
(Fig. 3c) invariant manifolds associated to one of the Lyapunov periodic orbits 
of the main family around Li and L2. Note that we plot the projection of the 
invariant manifolds on the configuration space (x, y). Our results show that the 
morphologies obtained depend on dynamical factors, that is, on the presence of 
homoclinic or heteroclinic orbits of the first intersection of the corresponding 
invariant manifolds. If heteroclinic orbits exist, then the ring of the galaxy 
is classified as ri?i (Fig. 3a). The inner branches of the invariant manifolds 
associated to 71 and 72 outline a nearly elliptical inner ring that encircles the 
bar. The outer branches of the same invariant manifolds form an outer ring 
whose principal axis is perpendicular to the bar major axis. If the model has 
neither heteroclinic, nor homoclinic orbits and only transit orbits are present, 
the barred galaxy will present two spiral arms emanating from the ends of 
the bar. The outer branches of the unstable invariant manifolds will spiral 
out from the ends of the bar and they extend azimuthally to more than 37r/2 
(Fig. 3d). If the outer branches of the unstable invariant manifolds intersect 
in configuration space with each other ^, then they form the characteristic 
shape of R2 rings (Fig. 3b). That is, the trajectories outline an outer ring 
whose major axis is parallel to the bar major axis. The last possibility is if 
only homoclinic orbits exist. In this case, the inner branches of the invariant 
manifolds form an inner ring, while the outer branches outline both types 
of outer rings, thus the barred galaxy presents an R1R2 ring morphology 
(Fig. 3c). 



* Note that they cannot intersect in phase space. 
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a - rRl ring structure b - rR2 rmg structure c - rRlRS ring structure d - Barred spiral structure 




-10 -5 5 10 -10 -5 5 10 -20 -10 10 20 -20 -10 10 20 



Fig. 3. Model rings and spiral arm structures. We plot characteristic examples of 
(a) rRl ring structure, (b) rR2 ring structure, (c) R1R2 ring structure and (d) 
Spiral arms. 

Wc also study the response of an axisymmotric component to a bar pertur- 
bation. We use the same axisymmetric potential and the same bar potential 
as in our models and the bar is introduced gradually, to avoid transients. 
Once the bar has reached its maximum amplitude, we consider a snapshot of 
the response simulation and we compare its morphology to the correspond- 
ing structure we obtain with our models. In Fig. 4 we show the results for 
the spiral arms case, by over-plotting the selected snapshot with our model. 
The white points represent the particle positions of the response study and 
the black lines are the unstable invariant manifolds. Note that the two match 
perfectly. 
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Fig. 4. Over-plot of the results obtained with the response simulation (white dots) 
and the invariant manifolds (black lines) in a model with spiral arms. 
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We compare our results with observational data (E. Athanassoula, M. 
Romero- Gomez, J.J. Masdcmont, C. Garci'a-Gomcz, in preparation) and we 
find good agreement. Regarding the photometry, the density profiles across 
radial cuts in rings and spiral arms agree with the ones obtained from obser- 
vations. The velocities along the ring also show that these are only a small 
perturbation of the circular velocity. 

Acknowledgements 

MRG acknowledges her fellowship "Becario MAE-AECI". 



References 

1. E. Athanassoula: Phys. Rep. 114, 319 (1984) 

2. E. Athanassoula: MNRAS 259, 328 (1992) 

3. B. Barbanis, L. Woltjer: ApJ 150, 461 (1967) 

4. R. Buta: ApJS 96, 39 (1995) 

5. G. Contopoulos, Th. Papayannopoulos: A& A, 92, 33 (1980) 

6. J.M.A. Danby: AJ 70, 501 (1965) 

7. W. Dehnen: AJ 119, 800 (2000) 

8. D.M. Elmegreen, B.G. Elmegreen: MNRAS 201, 1021 (1982) 

9. RB. Eskridge, J.A. Progel, R.W. Podge, et al.: AJ 119, 536 (2000) 

10. N.M. Ferrers: Q.J. Pure Appl. Math. 14, 1 (1877) 

11. D.E. Kaufmann, G. Contopoulos: A& A, 309, 381 (1996) 

12. G. Kuzmin: Astron. Zh. 33, 27 (1956) 

13. B. Lindblad: Stockholms Observatorium Ann. Vol. 22, No. 5 (1963) 

14. A. Lyapunov: Ann. Math. Studies 17 (1949) 

15. P.A. Patsis: MNRAS 369, L56 (2006) 

16. M. Romero-Gomez, J.J. Masdemont, E. Athanassoula, C. Garci'a-Gomez: 
A & A 453, 39 (2006) 

17. M. Romero-Gomez, E. Athanassoula, J.J. Masdemont, C. Garci'a-Gomez: 
A & A 472, 63 (2007) 

18. M.P. Schwarz: ApJ 247, 77 (1981) 

19. M.P. Schwarz: MNRAS 209, 93 (1984) 

20. M.P. Schwarz: MNRAS 212, 677 (1985) 

21. A. Toomre: ApJS 138, 385 (1963) 

22. A. Toomre: ApJ 158, 899 (1969) 

23. N. Voglis, I. Stavropoulos, C. Kalapotharakos: MNRAS 372, 901 (2006) 

24. N. Voglis, P. Tsoutsis, C. Efthymiopoulos: MNRAS 373, 280 (2006) 



